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Abstract. This article describes different approaches to investigate the behavior of parton evolution in 
the proton by exploiting various aspects of the hadronic final state produced in Deep Inelastic Scattering 
Events at HERA. 
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1 Introduction 

Measurements of the hadronic final state in deeply inelas- 
tic ep scattering (DIS) provide precision tests of quan- 
tum chromodynamics (QCD). At HERA data are col- 
lected over a wide range of the negative four-momentum- 
transfer Q 2 , the Bjorken variable x and the transverse mo- 
menta pt of hadronic final state objects. 

Studies of the hadronic final 
t state may be used to get insight 

into the dynamics of the parton 
cascade exchanged in low- a: lep- 
ton proton interactions. Fig. 1 
shows a generic DIS process in 
which a gluon from the proton 
undergoes a QCD cascade. The 
gluon interacts with the virtual 
photon via a hard photon-gluon 
process which can be calculated 
within perturbative QCD using 
an exact matrix element. The 
cascade itself represents an ap- 
proximation for an all order ma- 
trix element calculation and sev- 
eral prescriptions to describe the 
QCD dynamics within the cas- 
cade have been proposed. The 
most familiar one is given by 
the so called DGLAP evolution 
equations [1]. In these equations 
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Fig. 1. Diagram of a 

generic DIS process at 

low x. Here kt denotes 

the transverse momenta 

of the exchanged gluons, 

x g is the fractional mo- 
mentum of the parton 

taking part in the hard 

interaction and x is the 

Bjorken scaling variable 
a For the HI and ZEUS Collaborations 
b Talk presented at the EPS03, Aachen, July 2003 



the large logarithms in Q are 
resummed, neglecting log(l/x) 
terms. This practically corre- 
sponds to a strong ordering of 
the transverse momenta of the 



emitted partons, i.e. fc tj i << .. << k t , n « ■■ « Q 2 • 
DGLAP evolution is expected to break down at suffi- 
ciently low values of x, when the ordering no longer ap- 
plies. 

At very low values of x it is believed that the theoret- 
ically most appropriate description is given by the BFKL 
evolution equations [2]. These resum large logarithms of 
1/x up to all orders and impose no restriction on the or- 
dering of the transverse momenta within the parton cas- 
cade. Thus off-shell matrix elements have to be used to- 
gether with an unintegrated gluon distribution function, 
f(x,fA,kt), which depends on the gluon transverse mo- 
mentum kt as well as x and a hard scale fif. A promising 
approach to parton evolution at both low and large values 
of x is given by the CCFM [3] evolution equation, which, 
by means of angular-ordered parton emission, is equiva- 
lent to the BFKL ansatz for x — > 0, while reproducing the 
DGLAP equations at large x. 



2 Forward 7r°/Jet Cross Sections 

An extended parton ladder at low x leads to high k t par- 
tonic emission in the region close to the proton remnant 
('forward' region) to which measurements of jets and lead- 
ing particles, e.g. it are sensitive. Production of a forward 
7T° can be regarded as a refined version of forward jet pro- 
duction. In order to enhance the sensitivity to low- a; effects 
special selection cuts have been applied such as confining 
the ratio p\ , a j et )/Q 2 to values between 0.5 and 2 in- 
spired by a proposal in [4] . 



2.1 Forward tv° Cross Sections 

Inclusive forward tt° cross sections for transverse momenta 
Pt,7t° > 3.5 GeV are shown in Fig. 2 as a function of x 
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H1 Forward Jet Data 



HI preliminary 

DIR + RES, Q 2 +4p T = 
CCFM (CASCADE) 
D1R (LO DGLAP) 
2.0 < Q 2 < 8.0 GeV 2 



8.0 <Q <20.0 GeV~ 




Fig. 2. Forward 7r° cross section as a function of Bjorken-x in 
different regions of Q 2 compared with predictions of DGLAP 
and CCFM based QCD Models. 



for different regions of Q 2 . The data are compared with 
predictions of the Monte Carlo models RAPGAP [5] and 
CASCADE [6]. RAPGAP implements a QCD model based 
on Leading Order (0(a s ), LO) parton showers with ('DIR 
+RES') and without ('DIR') resolved photon structure. 
CASCADE is employed as an implementation of the CCFM 
evolution equation introduced above. The prediction by 
RAPGAP with a pointlike photon (DIR) is well below 
the data. A reasonable description of the data is achieved 
by including contributions from resolved virtual photons 
in the predictions and using a factorization scale of Q 2 + 
4p T 7T o. Note, that resolved contributions can be consid- 
ered to mimic a lack of ordering in transverse momen- 
tum as required for a DGLAP evolution scheme. CAS- 
CADE predictions based on the unintegrated gluon den- 
sity JS2001 [6] on the other hand undershoot the data for 
lower values of Q 2 . 
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Fig. 3. Forward jet cross section as a function of Bjorken-a; 
compared with NLO DGLAP QCD calculations and predic- 
tions by the CCFM Model CASCADE. 
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Fig. 4. Inclusive jet cross section as a function of Bjorken-a; 
compared with NLO DGLAP QCD predictions and DGLAP 
based QCD Models. 



3 Inclusive Jet Cross Sections 



2.2 Forward Jet Cross Sections 



Results complementary to the ones discussed in Sec. 2.1 
are obtained by studying jets in the same region of phase 
space. Jets are reconstructed with the longitudinally in- 
variant k t cluster algorithm [7]. Fig. 3 shows the forward 
jet cross section for transverse momenta pr. jet > 3.5 GeV 
as a function of x. The data are compared with NLO 
(0(a 2 )) QCD calculations performed with the program 
DISENT [9] and predictions by CASCADE based on two 
recent sets of unintegrated gluon distributions [8]. While 
results of the NLO QCD calculations are significantly be- 
low the data, the CASCADE prediction based on the set 
labelled J2003-1 is in good agreement with the data. The 
difference between the CASCADE predictions indicates 
the sensitivity of forward jet data to low- a; dynamics. 



In the following analysis of inclusive jet cross sections the 
restriction to the forward region and the kinematic con- 
finement introduced in Sec. 2 are removed, leading to a 
somewhat more general study of jet cross sections. In [10] 
it is outlined that the comparison of the measured jet cross 
sections with matrix element calculations including con- 
tributions up to 0(a s ) performed with DISENT lead to 
significant discrepancies between the data and the theoret- 
ical predictions. If, however, hadronic activity in the for- 
ward and the backward hemisphere is required and NLO 
(0(al)) predictions are employed a much better agree- 
ment between the data and the theoretical prediction is 
obtained, which is demonstrated in Fig. 4. The theoretical 
error represented by the hatched band is due to missing 
higher order contributions in the theoretical calculations. 
Predictions based on LO DGLAP parton showers, here 
represented by LEPTO [11], undershoot the data at low- 
er; while a good description of the data is obtained by the 
CDM [12] model as implemented in the event generator 
ARIADNE [13]. 
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4 Azimuthal Correlations between Jets 

Insight into low- a; dynamics can be gained from inclusive 
dijet data by studying the behavior of events with a small 
azimuthal separation, A(f>* , between the two hardest jets 
as measured in the hadronic center-of-mass system [14- 
16]. Partons entering the hard scattering process with 
negligible transverse momentum, k t , as assumed in the 
DGLAP formalism, lead at leading order to a back-to-back 
configuration of the two outgoing jets with A<j>* = 180°. 
Azimuthal jet separations different from 180° occur due to 
higher order QCD effects. However, in models which pre- 
dict a significant proportion of partons entering the hard 
process with large kt, the number of events with small 
A<j>* should also increase. 

Here we present a measurement of the ratio 



S 



J N diict (Ar,x,Q 2 )dAr 
jo 18 °° N diiet (Ar,x,Q 2 )dAr 



of the number of events N d ^ et with an azimuthal jet sep- 
aration of A(f>* < 120° relative to all dijet events. The 
observable was proposed in [16] and is considered to be 
directly sensitive to low-x effects. 

Fig. 5 presents the 5-distribution as a function of x for 
different Q 2 . It is of the order of 5% and increases with 
decreasing x. This rise of S is most prominent in the low- 
est Q 2 bin, where the smallest values of x are reached. 
The NLO dijet QCD calculations, resulting in an effective 
LO prediction for S, predict ^-values of only ~1% and 
show no rise towards low x. NLO 3-jet predictions by the 
program NLO JET [17] which lead to an effective NLO 
prediction for S give a good description of the data at 
large Q 2 and large x, but still fail to describe the increase 
towards low x, particularly in the lowest Q 2 range. Ac- 
cording to the discussion given in [18], predictions based 
on CCFM evolution as implemented in CASCADE lead to 
an improved overall agreement with the data in particular 
at low x and Q 2 . 
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